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Abstract. We estimate the distribution of the eigenvalues of a family of time-frequency localization operators 
whose eigenfunctions are the well-known Prolate Spheroidal Wave Functions from mathematical physics. These 
operators are fundamental to the theory of bandlimited functions and have applications in signal processing. 
Unlike previous approaches which rely on complicated formulas for the eigenvalues, our approach is simple: We 
build an orthonormal basis of modulated bump functions (known as wave packets in time-frequency analysis) 
which approximately diagonalizes the operator of interest. 


1. Introduction 


This paper concerns the distribution of eigenvalues of time-frequency localization operators (TFLOs) of the 
form 


( 1 ) 


Ti,jf = RiPjRif 


where I, J are compact intervals, and Ri : L^(R) —L^(R), Pj : L^(R) —L^(R) are associated projection 
operators in the time and frequency variables, i.e. 


( 2 ) 

( 3 ) 




f{x) X el 
0 X eR \ I, 


{Pjf){x) = [ [ /(2/)e-2--(^-)dydcc. 
J.j Jr 


The eigenfunctions of T/ j are the restrictions to I of Prolate Spheroidal Wave Functions (PSWFs), which 
arose first in the study of the Helmholtz equation in mathematical physics IITSl . Due to the cormection 
with TFLOs, the PSWFs provide an optimal basis for the representation of bandlimited fimctions on an 
interval, as observed in Landau-Pollack IHIIU, Slepian flOlITTl . and Slepian-PoIIack IIT2l . More recently, 
PSWFs have been applied to produce quadrature formulas ||8l and interpolation schemes ||9| for computing 
with bandlimited functions. 

As demonstrated in m, the eigenvalues of Tj j display a concentration phenomenon in the asymptotic 
limit as the parameter A := |/| • | J| tends to infinity: For e > 0, there are approximately A eigenvalues in the 
interval [1,1 — e], approximately ClogjA) • log(l/e) eigenvalues in the interval (e, 1 — e), and the remaining 
eigenvalues form a sequence tending to zero at an exponential rate. This estimate is asymptotic, meaning 
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that it is guaranteed only for some sufficiently large A; no quantitative bounds were known for any finite A 
until recently: A quantitative upper bound on the eigenvalue sequence was proven in Osipov |[7|. 

This paper's aim is to present an alternate method for estimating the eigenvalues of a TFLO using tech¬ 
niques from time-frequency analysis. We will prove quantitative upper and lower bounds by giving an 
example of a basis of time-frequency wave packets (the local-cosine basis) that approximately diagonal¬ 
izes the TFLO. Our results are sub-optimal by a single factor of log(l/e). We hope that the methods given 
here can be used to study localization operators associated to domains in higher-dimensions. We leave this 
investigation for a future work. 


1.1. Statement of main results. We fix intervals J = [—1/2, 1/2] and I = [—D/2, D/2\, where D >2. 

The operator T = T/, j defined in |(l]l is compact and positive-semidefinite, and its operator norm is at 
most 1. We write as the direct sum % 0 ker(r), where H = ker(T)-'*; clearly, L^(K \ /) is a subspace 

of ker(r). Thus, we may regard T as an operator on L^{I). The spectral theorem for compact operators 
implies that the spectrum of T is discrete and the positive eigenvalues of T form a non-increasing sequence 
{Afe} in (0,1], with A^ —0 as /c oo. Furthermore, there exists an orthonormal basis {'0fc}fceN for the 
subspace % = ker(T)'^ of L^{I) satisfying the eigenvalue equation for all /c S N. 

The position of the eigenvalue of T closest to A = 1/2 is described in the next result from ||4|. 


Theorem 1. We have 


( 4 ) 


A[_d]_i < 2 - ^[^ 1 ' 


Here, we write [a;] to denote the integer part of a real number x. 


Our main result bounds the number of eigenvalues contained in an interval centered about A = 1/2. 

Theorem 2. For each 77 G (0,1/2] there exists a constant A,, > 1 such that the following holds. 

Given e G (0,1/2) and D >2, define 

K = Ar,- (log (log(i?) log(i? ■ e-1). 

Then 

{fc G N : Afe G (e,l-e)} C [D-K,D + K]. 

Remark 1. In 0, a one-sided bound {fc G N : A^ > e} C [0,D + K'] is proven, with K' = Clog(Zl)^ log(l/e). 
Our result has an improved dependence on D, but sub-optimal dependence on e. The optimal bound is expected to be 
K" = C log(Zl) log(l/e), as predicted by the asymptotic analysis in |3|. 

Our strategy for proving Theorem |2| will be to construct an approximate eigenbasis for the operator T. We 
define a family of translated and modulated bump functions (the local-cosine basis). The result follows 
by counting the number of basis functions whose "time-frequency profile" is localized inside the rectangle 
/ X J in phase space. 
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The outline of the paper is as follows. In Section|2]we present the functional analysis lemma that will be 
used to control the eigenvalues of our operator. In Section |3] we construct a local-cosine basis and present 
its basic properties. We prove the key energy estimates for this basis in Section|4l We conclude the paper in 
Section|5]with the proof of Theorem|2l 


1.2. Acknowledgements. We are grateful fo Mark Tygert for useful discussions and for his comments dur¬ 
ing the preparation of fhis work. We are also grafeful fo Rachel Ward for proofreading an early version of 
fhis paper. 


2. The Main Lemma 

Let "H be a complex Hilbert space, and let {4>k}k&x be an orthonormal basis for "H. Denote the inner 
product on V. by (•, •) and the Hilbert norm on H by ||(p|| := {ip, ip) for tp e 

Let T : "H -A H be a compact, positive-semidefinite operator, with operator norm at most 1. By the 
spectral theorem for compact operators, there is a sequence 1 > Ai > A 2 > • • • > 0, and a complete 
orthonormal basis for fhe subspace %' = ker(T)-*“ of 'H, satisfying the eigenvalue equation T(p£ = 

X( ■ ipi for all € N. For e € (0,1/2), let = M^{T) denote the number of eigenvalues (counfed with 
multiplicity) that belong to the interval (e, 1 — e). 

Lemma 1. Let {(j)k}k^i be an orthonormal basis for LL. Assume that X is the disjoint union Iq U U I 2 , where 

(5) ^ ||T</)fe|p+ ^ <£3. 

/cGX2 

Then we have < 2 • #(Ii). 

Proof. Note that Mg is the dimension of the subspace Se '■= span{(p£ : A^ € (e, 1 — e)} of XL. 

Consider the orthogonal projection operator tt^ : S'g. Note that T and tt^ commute, since is 

spanned by a collection of eigenvectors of T. 

By definition, we have e\\(j)\\ < \\T(j)\\ and e||0|| < \\T(j} — fW for all f G S^. Hence, 
e- hefW < ||r7rg(/|| = ||7rgT(/|| < \\T(j)\\, and 

e • < WTiTef - = \\TTe{T(j)- /))|| < \\T f - (j)\\ for all feXL. 

Using f = fk (k G Xq U X 2 ) in the above estimates, we see that 

• hefkf < Y + Y ^ 

fcgioUi2 fceio fee2:2 

Therefore, 

(6) Y^ < e < 1/2. 

A:GXoUX2 
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We may assume that = dim(S'e) > 1, for otherwise the conclusion of the lemma is trivial. Then, by 
the Parseval identity and the fact that is an orthogonal projection operator for S^, we have 

llV'f = forall V' € Se. 

feel fcgi 

We fix an orthonormal basis for and sum the previous estimate over all ip belonging to that basis. Thus 
we obtain 

(7) dim(5',) = ^ ||7re<^fef, 

fcGl 

where here again we have used the Parseval identity to simplify the right-hand side. Recall that I is the 
disjoint union Iq U U X 2 . Since dim(S'e) > 1, we learn from ([61 and ([Tjl that 

(8) ^ > dimjS'e) - 1/2 > (1/2) • dim(S'e). 

feGli 

Finally note that ||7r£(/fe|p < ||(/fe||^ = 1 for any k G Xi. Thus ll8]l yields = dim(S'e) < 2 • #(Xi), as 
desired. □ 


3. Local Trigonometric Bases 

In this section we exhibit a smooth compactly supported cutoff function whose Fourier transform has 
near-exponential decay. We follow an approach found in 12. Using this cutoff function we construct an 
orthonormal basis for X^(/) (/ a compact interval) consisting of modulated bump functions. This is the 
local cosine basis of Coifman-Meyer jT]. 


3.1. A cutoff function. Fix an integer m > 1, and define 


(9) 


a{x) 


™g-(a:+l) ”* 

0 


X G (-1,1) 

X G {—00, 1] U [1, 00). 


Write f to denote the /c-fold derivative of a frmction / : R R. 


Lemma 2. We have 

\D'^a{x)\ < (16m)'= • fork>0, x€ R. 

Here our notation is that 0° = 1. 


Proof. First observe that D^{e ^ "^) (a; > 0) is equal to the sum of 2^ terms of the form 

where w is a real number, and j, r are integers satisfying j + r = k and |w| < ■ [{m + l)j -I- rf. This 

statement is easily proven by induction on k. 
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Using the estimate y^e ^ < R^, for y,R> 0, we obtain 


\Fw,j,rix)\ < ■ [(to + l)j + rY ■ [[(to + l)j + r] • to 

= • [[(to + l)j + r] . TO-i] (-+i)'h-+U-- 

<m’^ [(l + TO-i).fc]^'+'"''^'" 

< (4TO)'=-fc(i+'"”')-''. 


We conclude that |£)^(e ^ )| < (Sm )’^for a; > 0. Hence, by the Leibniz rule we have 


|U^’a(a;)| < 2^ • max 
0<k'<k 


sup 


sup [£)'=-'='(e-i/“")| 

xG(0,2] 


xe(o, 2 ] 


< 2 ^ max 

0<k'<k L 

< (16to)'“' • 


2^^'rrY'{k — 
for all a; S M. 


This completes the proof of the lemma. 


□ 


We define 

( 10 ) ^(a;) := • f a{y)dy. 

^}^ady J_^ 

Since a{y) is an even function (see l|9)), we have 

(11) A(x) + A{—x) = TT/2. 

Now let 9{x) := sin(A(a:)). From llTT]| we deduce that 0(—a:) = cos(A(a:)). Therefore, 

( 12 ) 0 ^(a;) + 0 ^(—a:) = sin^(A(a;)) + cos^(A(a:)) = 1 for any x € M. 

Since a{x) = 0 for x € K \ [—1,1], we have A{x) = 0 for x < —1, and A{x) = 7 r /2 for x > 1 . Hence, 

(13) 9{x) = 0 for X < —1; 9{x) = 1 for x > 1. 

The next lemma provides an estimate on the size of the derivatives of 9. 

Lemma 3. Let F ■. C ^ Cbe an entire function, and let f -.R ^Rbe C°°. 

Assume that there exist C >1 and 7 > 1 such that |H^/(x)| < ■ k^^ for all k > 0. 

Then there exists Cq >1 determined by C, 7, and F, such that \D^ [F{f{x))] | < Cq • k'^’^ for all k >0. 

A proof of Lemma |3] is given in 121 . 

From the definition (|9l it is clear that / a(y)dy > > 1/16. We apply Lemma |2] to bound the 

derivatives of A{x) defined in (HOl l. Thus we obtain 

|H'=yl(x)| < C'^ ■ for all fc > 0 . 
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We apply Lemma |3] to the functions F{z) = sin( 2 :) and f{x) = A{x). Thus, for each m > 1 there exists 
Cm > 1 such that 


(14) 


\D^0{x)\ < Ct ■ for all k>0. 


We summarize conclusions dT^ . dlSll . and ifMll , in fhe following result. 

Proposition 1. Given a real number 77 > 0 there exists a C°° function 0 : M ^ R satisfying 

(a) 9{x) = Ofor x < —1, and 6{x) = Ifor x>\. 

(b) 0^(x) + 9'^{—x) = Ifor all a; G R. 

(c) \D^6{x)\ < for all integers k >0 and all a; G R. Here, > I is a constant determined by tj. 

To obtain the previous result from (Il2]| , ifTSll , and ifMll . we choose m > r]~^. 

3.2. Whitney intervals. By a dyadic interval we mean an interval of the form {k ■ 2“^, (fc +1) • for k.l^'L. 

Let I = [—D/2, D/2] for D > 0. The Whitney decomposition of / is a collection TV = consisting of 

dyadic intervals. Its basic properties are as follows. 

(Wl): The intervals Ij {j G J) are pairwise-disjoint, and I = Ujg • 

(W2): We have \Ij \ < dist{lj,dl) < 5 ■ \Ij \ for all j G J. 

For a construction of the Whitney decomposition, see Ifl4l . 

3.3. The local cosine basis. Let TV = {Ij}jGj be the Whitney decomposition of an interval I. We write 

Ij = {xj,xj +Sj]ioTj G J. 

We choose positive real numbers pj and rj/ satisfying rjj + 77 ' < - Then define a C°° function 0^ : R 

R by the formula 



(15) 


Part (a) of Proposition[T]implies that 9j is supported on [xj — -^Sj, Xj -I- ^<^ 7 ] ■ Since dist(/ 7 , dl) > \Ij \ = dj 
(see (W2)), we have 


(16) 


supp 9j c I for all j G J'. 


Denote 


r = {O', k) ■. j ej, k ez, k> 0 }. 


For any (j, /c) G F we define ^{j,k) & by 


(17) 


^■(j, 7 c)(a;) = C( 7 ,fe) • 6j • 9j{x) cos (tt • 5 C . (fc + 1 / 2 ) • {x - xj)) , 


where 


(18) 
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This definition of C(j fc) ensures the normalization condition ||l 2 (/) = 1. 

Note that 9j > 1 on [xj + and that 6j is supported on [xj — -^Sj^Xj + Therefore, the 

value of the integral term in the parentheses in (flSt is between and 25j. We conclude that 

(19) i < too. 


A theorem of Coifman-Meyer ||l] states that {‘l’(j,fc)}(j,/c)er is an orthonormal basis for T^(/) for an ap¬ 
propriate choice of the constants rjj and r;' which satisfy 


( 20 ) 


1 

too 




This is often called the local cosine basis or the Coifman-Meyer basis. We note that the construction in 
ID uses a different cutoff function 9. However, the proof of orthonormality requires only the conditions 
foimd in parts (a) and (b) of Proposition [l5l Consequently, the arguments in m establish orthonormality 
for the basis constructed here. We note that our cutoff function 9 satisfies the derivative bounds in part (c) 
of Proposition[T5l which will be used later. 


We henceforth assume that rjj and 77 ' satisfy ll20ll and are chosen so that {^{j,k)}{j,k)Gr is an orthonormal 
basis for 

We write / or T{f ) to denote the Fourier transform of a fimction / G L^(R), defined via the formula 

m = [ /(x)e-2"“«dx. 

Jr 


We require the following lemma from 121 . 


Lemma 4. Let 9 G (^““(R). Let C > 1 and 6 >lbe such that 

(a) 9 is supported on [—1,1], and 

(b) \D'^9{x)\ < for all k >0 and a: G R. 

Then |0(^)| < Aexp(—a • for all ^ G R, where a, A > 0 depend only on C and 5. 


We define ipj{x) := 9j{x ■ Sj +Xj) for j G J, which can be rewritten as fj^x) = 9 (^x ■ ■ 9 y{l — x) ■ 

(see GUI). From part (c) of Proposition [T] and since G we conclude that \D^'ijjj{x)\ < 

for k > 0 and a; G R. By applying Lemma |4] we learn that 

IfjiOl < ■ exp (-a,) • 

Because of the scaling relationship between ipj and 9j and simple properties of the Fourier transform, as 
well as the bound 1 — 77 < (1 -I- 77 )“^ (77 > 0), we conclude that 

(21) 1 ^ 7(01 ^ ' ^3 ' exp(—a,, • for w G R. 
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Using the formula iflTll and the scaling/translation/modulation properties of the Fourier transform T, 
we have 



^ • exp (—7ri • (fc + 1/2) • ^^ ^ • Xj) 
I ■ exp {iri ■ {k + 1/2) • S~^ ■ Xj) . 





In particular, thanks to (|2T]| and we have 



for constants a,, > 0 and Crj > 0. 


For k G Z>o and j & J we denote 


— (2fc + 1) • (4(5j) 


If we let Bn{^) := exp (—a,, • ''), then the bound (l22t states that 


(23) 


|^($0-.fc))(e)l < 5 ]^^ ■ [Br, (S, • (e - ^Jk)) + (s, ■ (e + e,fc))]. 


4. Energy estimates 


Recall that J = [—1/2,1/2] is the frequency localization interval and I = [—D/2,D/2] is the time local¬ 
ization interval. We decompose I into its Whitney decomposition W = We write Ij = {xj , Xj -i 6j]. 

In the previous section we defined an orthonormal basis {^{j,k)}j&j,kei,>o Recall that ‘Fq is 

supported on the interval I* = {xj — -^Sj,Xj -i Moreover, the Fourier transform of is (nearly) 

exponentially concentrated about the frequencies ^ = i^jk in the sense of the bound ll23ll . In this section we 
will derive the main energy estimates on our basis. 

Let s > 1 and (5min S (0,1) be parameters, which will be determined in the next section. 

We write X = 0{Y) to indicate the inequality |X| < C Y, where C is a constant independent of all 
parameters. We write X = 0{Y) to indicate the inequality |X| < Crj -Y, where C is a constant depending 
only on the parameter t]. 

Consider the basis functions ^{j,k) {k G Z>o) associated to a fixed Whitney interval Ij G W. We partition 
this collection into three groups by partitioning the index set Z>o as follows: 


(24) 

(25) 

(26) 


:= {k G Z>o : dist((2fc -E l)/(4dj), R \ J) > s • J-^}, 
_ 1 ^ g . dist((2fc + l)/{4Sj),dJ) < s ■ d-i}, 
:= {k G Z>o : dist((2fc + 1)/{U,), J) > s ■ Sj^}. 


Lemma 5. We have 


. #(£)-)= d, +0(s) ffd, >s, 
. #(£)-) = 0 ifSj < s, 
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• < lOs. 


Proof. The numbers ^jk = (2fc + l)/(4(5j) {k € N) form an evenly-spaced grid of width ^ starting at 
^jo = Recall that J = [—1/2,1 /2]. Clearly there are at most Sy +1 many indices k € N satisfying G J. 
Furthermore, if fc < Sj — 2s — ll2 then < 1 _ and thus dist(^jfc,R \ J) > s ■ S~^. Therefore, 

we have 

Sy -2s- 1/2 < #(£*“'") < Sy + 1. 

This implies the first bullet point. 

If Sy < s then s • > 1. There are no points whose distance to K \ [— 1 /2,1/2] is greater than 1. Thus, in 

this case = 0. 

The spacing between consecutive numbers ^jk is equal to 2 ^. Thus, at most 2s -|- 1 of the ^jk lie in an 
interval of width s • about the boundary point 1/2 G dJ. The same is true for the boundary point 
— 1/2 G dJ. Therefore, < 2 ■ (2s -I- 1) < 10s. 

□ 


We partition the index set T = x Z>o into three components: 

(27) Tiew = {(j, k) : Sy > S„,in, k G 

Tmed = {(j, k) : Sy > Smin, k G 

Thigh = {(/, k) : Sy > dniin, k G U {(j, k) : Sy < S^in, k G Z>o}. 

Lemma 6. For a numerical constant C >0,ive have #{Tmed) < Cs ■ log(i7/(5min)- 


Proof. Lemma |5] implies that 

#(rn.ed)= ^ *{£T^)< Y. 10®- 

Property (W2) implies that the number of Whitney intervals Ij for which 5j > Jmin is boimded by 
C'log(diam(/)/(5min) = Clog(£>/5inin)- Hence, 

#(rmed) < Cs ■ log(Zl/(lmin), 

as desired. 


□ 


We will now prove that Fourier transform of a basis function indexed by {j, k) G Flow or (j, k) G Thigh 

is sharply concentrated on J or M \ J, respectively. 
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Lemma 7. There exist constants C,c> 0 determined by y such that 


(28) ^ ll■^(^O^fe))lli=(E\J) < C'exp(-c • s^-^) ■ \og{D/S^,^) 

(29) ^ < Cexp(-c • s^-^) ■ \og{D/5n,in) + CS^in 


Proof. Recall our notation: ^jk = for (j, fe) G F = J' x Z>o. 

For j G J with Sj > (5min/ we define 

(30) 47 := {k G Z>o : dist(ejfe,R \ J) G [s • 2^/Sj,s • 2^+^/dj)}, and 

(31) :={fcGZ>o :dist(0fc,d) G [s-2V<5j,s-2^+V(5j)} for£GZ>o. 

Note that 4°" = Uf>o 4°J ^ ^£>0 lEU-lESll. 

The spacing between and for distinct k, k' G Z>o is at least Thus, a counting argument shows 


that 


(32) 

#(>C“f) < 10s • 2^ 

and 


(33) 

#(£‘°J) < 10s • 2^. 


From (l23ll . for any k G Z>o we have 

ll-^('J’a,fe))lli^(R\J) <c [ [Br,{S, ■ (e - e,fc)) + Br,{6j ■ (e + e,7)] 'rfC- 

Jr\j 

Since IR \ J = (—oo,l) U (l,c)o) is symmetric about the origin, we deduce that 

ll•^(<^>0^fe))lli2(R\J) <C [ <5, • B,{S, • (e - 

dR\J 

<c [ Br,{efde 

where the second inequality relies on the change of variable 7 = — note that G R\ J => |^ > 

Sj ■ dist(^jfc, R \ J). Thus, by the definition Bn{^) = An ■ exp(—a,, • we conclude that 

(34) ll•^(^(i.fc))llL2(R\J) < Cexp(-c- [,5j ■ dist(^jfc, R \ J)]^”’*) 

for constants c,C > 0 that depend only on y. For k G we have dist(^jfc,R \ J) ~ s • 2^fSj, where we 
write A ^ i? to mean that cA < B < CA for some constants c, C. Thus, (l34ll implies that 

(35) ll■^(‘&0.fe))lli=(R\J) < Cexp (-C- [s • 2']'“'') for all k G 
The method used to prove (l34ll also shows that 

(36) ll•^(^ 0 .fe))lli 2 (J) < Cexp^-c- •dist(Cjfe,d)]^“'') 
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for constants c,C > 0 depending only on 77 . For k G £’^f' we have dist(^jfe, J) ^ s ■ 2^/Sj. Thus, 
that 

(37) 


implies 


ll•^(^ 0 ■.fc))lli 2 (J) < C'exp^^-c- [s-2^]^ for all fc 


G . 


We write = U^>o ll33l and (l35l . we learn that 

CXD 

(38) E E = E E E \\^i^u^))\\lH^\j) 

fQ^£}ow j'Sj^Sjnin ^—0 

00 

< E El'^'®2^'‘3'exp(-c-[s2V“’ 

j:5j >l5min (=0 

< E C • exp( —C • 5^“’') 

j'Sj ^l5niin 

< C exp(-c • s^”’') log(diam(/)/(5i„in). 


In view of the definition of Plow in ll27ll . this completes the proof of I 

Next we prove ll29t . We write = U^>o From (I32t and ll^ we have 

00 

(39) E E \\n^u,>^))\\Uj)= E E E 

00 

i N E 10s2^ • C'exp(—c • [s2^Y 
j(5j>5min 

< E C ■ exp(—c • s^”’') 

3'^j 

< C exp(-c • s^”’') log(diam(/)/(5i„in). 


Alternatively, suppose that j G is such that 5j < ^min- Then (l23t implies that 

(40) E \\^ i ^ im )\\ hij ) < T.hr • ^^-^ 3k )) + B ^{ 5 , ■ (e + e,fe))]'de. 

fc>0 k>0'''^ 

Switching the order of summation and integration in (l40t and using the fact that the interval J = [—1/2,1/2] 
is symmetric about the origin, we have 

( 41 ) E \\H‘^U,k))\\hiJ) <c jY.^r - ^3k))]'di. 

k>0 4jj.>g 

Since Sr;(0 is smooth, bounded, and rapidly decaying as ^ 00 , we can compare a Riemann sum with an 

integral to prove the estimate 


A:>0 


Br^{z)‘^dz < C uniformly for all ^ G M. 
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Therefore, from drill we have 

fc >0 

By summing over all j G J with Sj < ^min/ we conclude that 

(42) ^ Y.\\^^^im\\Uj)<C J2 ^J<CSn.in, 

where the last estimate is a consequence of the Whitney conditions (Wl) and (W2) (see Section|32). 

In view of the definition of Thigh in (l27l l. we see that (l3^ and (|4^ imply the estimate ll29ll . finishing the 
proof of the lemma. □ 


5. Proof of Theorem[2] 


In the previous section we defined an orthonormal basis for depending on a pa¬ 

rameter r] O (0,1/2]. 

Let e G (0,1/2). Let s > 1 and dmin G (0,1) be parameters. We will choose s and Jmin in the following 
paragraphs. 

In the previous section we defined a partition of T as Tiow U Tmed U Thigh in terms of the parameters s and 
<imin; see (I24l l- (l26t and (1271 . 

For all / G PlanchareTs theorem implies that 

WTfWl.^i) = WRiPjRifWh^i) = WRiPjfWlm) 

^ ll^.//lli2(R) 

Similarly, 


11/ - T’/lli^(7) = 11/ - RiPjRifWh^i) = \\Rif - RiPjfWh^i) 

< Ui-Pj)f\\hm 



Thus 

(43) y] ii^4>(j^fe)iii2(/) + y wp^u.k) -^{j,k)\\L'^(i) 

< y m<^u,k))\\Uj)+ y ii-^(4^LT-))iii=(E\j) 

(J 5 ^^ ^high (J 5 ^£ ^low 

Lemma|3 

^ Cf) ■ 6Xp( C-q ‘ S ■ log(Z)/^min) “t“ Cq ' ^min; 


where Cq,Cq > Q are constants determined only by rj. 
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We are ready, at last, to state our assumptions on s and Jmin- We take 

^3 


^min ■ — 


(44) 


2Cr, 


and 


■ exp(-c^ ■ '') • log(£>/5inin) < 2 ^^ 


The second estimate is equivalent to 


Elementary algebra shows that it suffices to take 

(45) s:= A, - (log (log(Z?) 

for a constant determined only by rj. Now using ll44ll in ll43]l we see that 


Recall that T is equal to the disjoint union Flow U Fmed U Thigh- Thus, according to Lemma [TJ if we let At 
{k € N) denote the positive eigenvalues of T (arranged in non-increasing order), and if we define to be 
the number of eigenvalues of T in the interval (e, 1 — e), then we have 

Lenrunaf^ 

M, < 2 • #(r,„,d) < Cs ■ log(-D/^„,i„) 

= CAr, ■ (log (log(-D) • . log(2C,,-D • e'^) 

< a; • (log (log(i^) • • log(-D . e-i), 

for a constant A,, determined only by rj. We apply ||4} and conclude that {fc : S (e, 1 “ e)} C [D — 2M^, D + 

2Mf]. According to our upper bound on M^, this yields 

{k : Afe G (e, 1 - e)} C[D-K,D + K], 


where 

A = A" . (log (log(-D) . £- 1 ))'+"" ■ log(-D • £- 1 ) 
for a constant A" determined only by 77 . This completes the proof of Theorem|2l 


□ 


References 

[1] R.R. Coifman and Y. Meyer. "Analyse de Fourier adaptee a une partition par des intervals de Whitney," Colloq. Math. 63: 111-117, 
1992. 

[2] J. Dziubahski and E. Hernandez. "Band-limited wavelets with subexponential decay," Canad. Math. Bull. 41: 398-403, 1998. 

[3] H.J. Landau. "Eigenvalue Distribution of Time and Frequency Limiting," Journal of Math. Analysis and Applications 77: 469-481. 

[4] H.J. Landau. "On the Density of Phase-Space Expansions," IEEE Trans, on Information Theory 39 (4): 1152-1156 (1993) 

[5] H. J. Landau and H. O. Poliak, "Prolate spheroidal wave functions, Fourier analysis, and uncertainty - II," Bell System Technical 
Journal 40 (1): 65-84,1961. 

[6] H. J. Landau and H. O. Poliak, "Prolate spheroidal wave fimctions, Fourier analysis, and uncertainty - ni: the dimension of the 
space of essentially time- and band-limited signals," Bell System Technical Journal 41 (4): 1295-1336, 1962. 


13 



[7] A. Osipov, "Certain upper bounds on the eigenvalues associated with prolate spheroidal wave functions/' Applied and Compu¬ 
tational Harmonic Analysis 35 (2): 309-340, 2013. 

[8] A. Osipov and V. Rokhlin, "On the evaluation of prolate spheroidal wave functions and associated quadrature rules," Applied 
and Computational Harmonic Analysis, http://dx.doi.Org/10.1016/j.acha.2013.04.002,2013. 

[9] A. Osipov and V. Rokhlin, "Detailed analysis of prolate quadratures and interpolation formulas," preprint (available online). 

[10] D. Slepian, "Prolate sphereoidal wave functions, Fourier analysis, and uncertainty - IV: extensions to many dimensions; general¬ 
ized prolate spheroidal wave functions," Bell System Technical Journal 43 (6): 3009-3057,1964. 

[11] D. Slepian, "Prolate spheroidal wave functions, Fourier analysis, and uncertainty - V: the discrete case," Bell System Technical 
Journal 57 (5): 1371-1430,1978. 

[12] D. Slepian and H. O. Poliak, "Prolate spheroidal wave fimctions, Fourier analysis, and uncertainty - I," Bell System Technical 
Journal 40 (1): 43-63,1961. 

[13] J.A. Stratton, "Spheroidal functions". Proceedings of the National Academy of Sciences 21: 51-56 (1935). 

[14] E.M. Stein, "Singular Integrals and Differentiability Properties of Fimctions," Princeton U. Press, 1970. 


14 


